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Abstract —This paper deals with the problem of relay selec¬ 
tion in wireless powered cooperative networks, where spatially 
random relays are equipped with energy storage devices e.g., 
batteries. In contrast to conventional techniques and in order to 
reduce complexity, the relay nodes can either harvest energy from 
the source sigual (lu case of uncharged battery) or attempt to 
decode and forward it (in case of charged battery). Several relay 
selection schemes that correspond to different state information 
requirements and implementation complexities are proposed. The 
charging/discharging behavior of the battery is modeled as a two- 
state Markov chain and analytical expressions for the steady- 
state distribution and the outage probability performance are 
derived for each relay selection scheme. We prove that energy 
storage significantly affects the performance of the system and 
results in a zeroth diversity gain at high signal-to-noise ratios; 
the convergence floors depend on the steady-state distribution of 
the battery and are derived in closed-form by using appropriate 
approximations. The proposed relay selection schemes are gener¬ 
alized to a large-scale network with multiple access points (APs), 
where relays assist the closest AP and suffer from multi-user 
interference. 

Index Terms —Cooperative networks, relay selection, SWIPT, 
RF harvesting, stochastic geometry, outage probability. 

I. Introduction 

T he roll-out of the Internet of Things will lead to the 
massive deployment of sensor nodes and a vast amount of 
information exchange, making it impractical, even impossible, 
to individually recharge/control these devices on a regular 
basis. Wireless powered communication (WPC) is a promising 
energy solution for the future highly dense and heterogeneous 
networks with major impact on many different applications. 
It refers to communication networks where particular nodes 
power their operations by the received electromagnetic ra¬ 
diation. The fundamental block for the implementation of 
this technology is the rectifying-antenna (rectenna) which is 
a diode-based circuit that converts the radio-frequency (RF) 
signals to direct-current voltage 0, a. 

WPC appears in the literature in three basic architectures 
that refer to different application scenarios; a rigorous survey 
of WPC is given in a, a. The first main architecture is 
the wireless power transfer (WPT), where a dedicated RF 
transmitter wirelessly transfers power to devices a, a. In 
contrast to ambient RF harvesting, WPT can be continuous 
and fully controlled and therefore is attractive for applications 
with strict quality-of-service constraints. The second network 
architecture is the wireless powered communication network 
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(WPCN), where a dedicated RF transmitter broadcasts energy 
at the downlink and WPT-based devices transmit information 
at the uplink Q, 0. The third fundamental architecture 
is the simultaneous wireless information and power transfer 
(SWIPT), where the RF transmitter simultaneously conveys 
data and energy at the downlink devices a, noi. Due to 
practical constraints, SWIPT cannot be performed from the 
same signal without losses and practical implementations split 
the received signal in two parts, where one part is used 
for information transfer and another part is used for power 
transfer. This signal split can be performed in CD i) time 
domain i.e., time switching (TS), ii) power domain i.e., power 
splitting (PS), or iii) spatial domain i.e., antenna switching 
(AS). 

A special case of SWIPT with particular interest for the 
cooperative networks is the SWIPT with energy/information 
relaying. In this network structure, a batteryless relay node 
extracts both information and energy from the source signal, 
and then uses the harvested energy to forward the source signal 
to a destination. In CD, the authors study the performance of 
a three-node Amplify-and-Forward (AF) relay channel, where 
the relay node employs TS/PS to power the relaying link. This 
work is extended in ifTSll for a Decode-and-Forward (DF) relay 
channel and the throughput performance is analyzed in closed 
form for both TS/PS techniques. A three-node relay channel 
with direct link, which combines TS-SWIPT with the dynamic 
Decode-and-Forward (DDF) protocol, is analyzed in d.The 
work in ifTSl studies a multiple-input multiple-output (MIMO) 
relay channel, where the relay node employs AS technique for 
SWIPT; the optimal AS is formulated as a binary knapsack 
problem and low-complexity solutions are proposed. Recent 
studies re-examine the fundamental three-node SWIPT relay 
channel by assuming that the relay node has full-duplex (FD) 
radio capabilities m, C3. In ifTfill . the authors study the 
optimal beamforming for a FD MIMO relay channel where a 
relay node uses TS/PS in order to power the relaying link; an 
appropriate beamforming vector protects the relay input from 
the loop-channel effects. The integration of the loop-channel 
interference into the WPT process is discussed in 03, and 
the optimal power allocation and beamforming design at the 
relay node are derived. 

In contrast to the above studies which focus on the basic 
three-node relay channel, more complicated network structures 
appear in the literature. A network topology where multiple 
source-destination pairs communicate through a common DF 
relay node that employs PS-SWIPT is discussed in ifTSll . In 
HD, the authors study a two-way relay channel where multiple 
communication pairs exchange information through a shared 
PS-SWIPT MIMO relay node; the optimal beamforming vec- 
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tor that satisfies some well-defined fairness criteria is investi¬ 
gated. The work in 1^ introduces multiple SWIFT relays 
and studies the interference relay channel, when multiple 
sources communicate with their destinations through dedicated 
SWIFT relays. On the other hand, the authors in ETl study 
the relay selection problem when an access point (AF) com¬ 
municates with the destination through multiple FS-SWIFT 
relay nodes, which are randomly located in the space. This 
work analyzes the outage probability performance for three 
main relay selection techniques; i) random relay selection, 
ii) relay selection based on the Euclidean distance, and iii) 
distributed (cooperative) beamforming. The consideration of 
the spatial randomness in WFC studies is of vital importance, 
since RF harvesting highly depends on the path-loss e.g., ED- 
ll^ use stochastic geometry tools in order to study different 
WFC networks. 

A SWIFT-based relay extracts both information and energy 
from the source transmission. For single-antenna relays, the 
TS technique requires perfect time synchronization and suffers 
from interrupted information transmission, since dedicated 
time slots are used for power transfer. On the other hand, 
the FS technique dynamically splits the received signal at 
the relay node and requires appropriate circuits that increase 
the implementation complexity. The implementation of FS/TS- 
SWIFT relay channel is an open research problem and several 
practical issues should be resolved. 

On the other hand, the integration of an energy storage 
device (e.g., battery, capacitor etc.) at the relay nodes, which 
is charged by the received RF radiation, introduces another 
WFC-based relay structure with new potentials and challenges; 
this network architecture has not been fully explored. The 
work in ll24l introduces a three-node relay channel where the 
relay node is equipped with a multi-level battery that can be 
charged by the source signal; the optimal switching policy 
between WFT and information decoding at the relay node is 
investigated. This work is extended in ll^ for a scenario with 
multiple relay nodes where a single-best relay node is selected 
for relaying; however, this work does not take into account 
spatial randomness, which is critical for WFC networks. A 
WFC network with spatial randomness and batteries is studied 
in ll 2 ^ ; specifically, the authors study the performance of a 
large-scale cognitive radio network, where cognitive nodes 
charge their batteries from the primary transmissions. 

In this paper, we focus on a WFC cooperative network and 
we study the problem of relay selection, when relay nodes 
are equipped with batteries and are randomly located in the 
coverage area of an AF. We assume that the relay nodes hold 
two-state storage devices (e.g., batteries, capacitors) which 
are charged by the source transmission; relay nodes become 
active and can assist source transmission, only when are fully 
charged. Based on this setup, we re-design the relay selection 
schemes proposed in (for conventional batteryless relays) 
and we study their performance. We investigate several relay 
selection policies that have different channel state information 
(CSI) requirements and correspond to different complexities. 
Analytical results for the outage probability performance of the 
proposed schemes as well as simplified asymptotic expressions 
for the high signal-to-noise ratio (SNR) regime are derived. We 


prove that in contrast to the conventional case ED, where a 
diversity gain equal to one is ensured, the battery model results 
in a zeroth diversity gain at high SNRs. The associated outage 
probability floor becomes lower when the battery status is 
taken into account with regards to the relay selection decision. 
Furthermore, the relay selection schemes are generalized to a 
large-scale network with randomly located AFs, where multi¬ 
user interference affects relay decoding. To the best of the 
authors’ knowledge, the problem of relay selection for WFC 
cooperative networks with batteries at the relay nodes, has not 
been reported in the literature. 

Paper organization: We introduce the system model and the 
main assumptions in Section |II] In Section [HIl we present 
the relay selection policies and we analyze their outage 
probability performance. Section |IV] generalizes the relay 
selection schemes for a multi-cell network. Simulation results 
are presented in Section |Vl followed by our conclusions in 
Section |Vl] 

Notation: denotes the d-dimensional Euclidean space, 

\z\ denotes the magnitude of a complex variable z, P(X) 
denotes the probability of the event X, E(-) represents the 
expectation operator; 7 ( 0 , x) denotes the lower incomplete 
gamma function ll26l Eq. 8.350], and 2 Fi{o.,b;c-,x) is the 
Gaussian or ordinary hypergeometric function ll26] Eq. 9.100]. 

II. System model 

We consider a single-cell where an AF communicates with 
a destination, D, via the help of a set of relay nodes Ri, 
whose locations are determined from a Foisson point process 
(FFF). All nodes are equipped with a single antenna and the 
relay nodes have WFT capabilities. The AF is backlogged and 
transmits with a fixed power P and a spectral efficiency tq 
bits per channel use (BFCU); the AF’s transmitted signal is 
the only WFT source for the relay nodes. Fig. [D schematically 
presents the system model. 

A. Topology 

We consider a disc, denoted by V, where the AF is located 
at the origin of the disc and the radius of the disc is p. The 
relay nodes are located inside the disc and their locations form 
a homogeneous FFF <I> (inside V) of density A; N denotes the 
number of the relays. The distance between the AF and the 
destination is denoted by c?o and there is not a direct link 
between the AF and the destination e.g., due to obstacles or 
severe shadowing. We assume that di denotes the Euclidean 
distance between AF and the i-th relay Ri with 1 < i < N, 
while Ci denotes the distance between Ri and D. 

B. Cooperative protocol 

The relay nodes have half-duplex capabilities and employ 
a DF policy. Time is slotted and communication is performed 
in two orthogonal time slots. Specifically, in the first time 
slot i.e., broadcast phase, the AF broadcasts the signal to 
the relay nodes. In the second time slot, a single relay (or a 
group of relays), which successfully decoded the source signal, 
forwards the signal to the destination according to the rules of 
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Fig. 1. The network model; an AP communicates with the destination D 
via the help of relays. 


the considered DF-based relay selection scheme. It is worth 
noting that the selected relay transmits at a rate tq (same with 
the rate of the source). 

C. Battery model 

Each relay node is equipped with a single energy storage 
device (e.g., battery, capacitor, etc.), which can store energy 
for future use. At the beginning of the broadcast phase, the 
battery can be either fully charged or empty (two states) ll^ 
Sec. II.A]; a relay with a charged battery is active and can par¬ 
ticipate in the relaying operation, while a relay with an empty 
battery is in harvesting mode and uses the AP’s signal for WPT 
purposes. We assume that the battery can be used either for 
charging or transceiver operations (a single battery/capacitor 
cannot be charged and discharged simultaneously lIZTl l. Fig. 
| 2 ] schematically shows the battery model and the associated 
switching mechanism. 

The WPT process is based on the AP’s transmissions and 
WPT from relaying signals is negligible i.e., relay nodes 
transmit in a much lower power than the AP and thus 
cannot satisfy the rectennas’ sensitivity requirements. Since 
the battery has only two discrete states (full or empty), an 
empty battery can be fully charged when the input power 
is larger than the size of the battery. Let = T'P be 
the size/capacity of the battery with T* < 1. The battery 
of the selected relay (if it is charged) is connected to the 
transceiver’s circuit and is fully discharged at the end of 
the second/relaying time slofl If the decoding is successful, 
most of the available energy is used for transmission i.e., 
the relay node transmits with a fixed power Pr (transmission 
dominates the total energy consumption). If the decoding is not 
successful, the available energy is mainly consumed in order to 
operate the basic receiver’s components (e.g., RF electronics, 
signal processing, etc.), transmit a Negative-Acknowledgment 
signal to indicate unsuccessful decoding as well as for static 
maintenance operations (e.g., cooling system). 

Markov chain model: The charging/discharging behavior of 
the battery can be represented by a finite-state Markov chain 
(MC) with two states {sq, si}, where the state sq indicates that 
the battery is empty and the state si that the battery is charged. 

'This assumption ensures that the state-transition probability matrix is 
independent of the channel conditions (path-loss, fading) and simplifies the 
analysis. 



Fig. 2. The battery model; the battery is either connected to the WPT device 
for charging or to the transceiver device for discharging. 


For the sake of presentation, we introduce the function S{Ri), 
which returns the state of the battery for the i-th relay i.e., 
S{Ri) = Sj means that the battery of the z-th relay is in 
the Sj state, where j € {0,1}. The state-transition probability 
matrix 11 can be written as 


n = 


1 — TTo TTo 
TTl 1 - TTl 


( 1 ) 


where ttq denotes the probability that the input power is greater 
than the battery size of the relay, and tti is the probability that a 
charged relay node is selected for relaying. Both probabilities 
are defined in the following section for the proposed relay 
selection schemes. If 77 = [770 pi] denotes the stationary steady- 
state probability vector of the MC, by using basic queueing 
theory, we have 


7711 = 77 . 


( 2 ) 


The solution of the above system of linear equations gives the 
steady-state distribution of the battery, which is equal to 

TTl TTo 


Vo = 


Vl = 


(3) 


TTo -h TTl TTq -h TTl 

Extension to a multi-state battery model: If a relay node be¬ 
comes active only when is fully charged, a battery model with 
L-l -2 states follows the same analysis with the two-state case; 
it only affects the computation of the steady-state distribution. 
Specifically, the battery of size is discretized in L-l-2 energy 
levels Ci = iPr/{L 1) with z = 0,..., L -|- 1. In this case, 
the harvesting threshold is smaller than the size/capacity of the 
battery and is equal to the difference between two successive 
energy levels i.e., Pr/{L-\- 1). We define L-l-2 corresponding 
energy states, Si, with z = 0,1,..., L-l-1 and thus the battery is 
in state Si when its stored energy is equal to a. If pij denotes 
the transition probability P{si —>■ Sj}, the state-transition 
probability matrix 11 has a dimension (L -|- 2) x (L -|- 2). A 
steady state 77 exists and can be calculated by solving the linear 
system in (|2]i. We note that the probabilities pij with i ^ L-|-2 
depend on the channel statistics and can be calculated using 
the procedure in ll24l . On the other hand, PL+ 2 ,j = 0 for 
0 < j < L -f 2, while the transition probabilities pL+ 2 , 0 , 
Pl+ 2 ,l +2 depend on the relay selection schemes and follow 
the same analytical steps with the two-state case. 


D. Channel model 

We assume that wireless links suffer from both small-scale 
block fading and large-scale path-loss effects. The fading is 
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Rayleigh distributed so the power of the channel fading is an 
exponential random variable with unit variance. We denote by 
hi and gi, the channel coefficients for the links between AP 
and Ri, Ri and D, respectively. The path-loss model assumes 
that the received power is proportional to 1/(1 + (i“) where 
d is the Euclidean distance between the transmitter and the 
receiver, a > 2 denotes the path-loss exponent and we define 
5 = 21 a. The considered path-loss model ensures that the 
path-loss is always larger than one for any distance i.e., 1 -|- 
(i“ > 1, even if d < 1. The instantaneous fading channels 
are known only at the receivers, except if otherwise defined. 
In addition, all wireless links exhibit additive white Gaussian 
noise (AWGN) with variance tr^; m denotes the AWGN at 
the z-th node. 

The proposed relay selection schemes are analyzed in terms 
of outage probability i.e., the probability that the destination 
cannot support the target spectral efficiency. If C{x) = 
^ log 2 (l + a;) denotes the instantaneous capacity for a wireless 
link (one-hop transmission) with SNR x, the associated outage 
probability is given by P{C'(a;) < rp} = P{a; < e} where 
e = 2 ^’’“ — 1; to simplify the notation, we define 5 = ecr^/P. 

III. Relay selection and pereormance analysis 

In this section, we present several relay selection schemes 
that correspond to different complexities. The outage probabil¬ 
ity performance as well as the diversity gain of the proposed 
schemes are derived in closed form. 


A. Random relay selection 


The random relay selection (RRS) scheme does not require 
any feedback about the battery status or the location of the 
relay nodes and selects a relay in a random way. It corresponds 
to a low implementation complexity and is appropriate for 
networks with strict power/bandwidth constraints. Without loss 
of generality, we consider that the z-th relay is selected to 
assist the source. If the z-th relay is fully charged, it attempts 
to decode the source signal and acts as a relay in case of 
successful detection. If the z-th relay has an empty battery, it 
switches to the harvesting mode and uses the received signal 
for WPT purposes; in this case, the relay remains inactive 
during the cooperative slot and an outage event occurs. On the 
other hand, the non-selected relays with empty batteries switch 
to harvesting mode and use the source signal for potential 
charging. During the first time slot, the received signal at the 
z-th relay can be written as 

y* = + (4) 

where s denotes the source signal with normalized power. If 
the z-th relay is active, the associated received SNR is equal 
to 


SNR, 


(i + d/V 


(5) 


If the z-th relay is inactive (empty battery), the input power at 
the WPT device is equal to 


Ph = 


PM^ 

(i + df)’ 


( 6 ) 


where C, denotes the WPT conversion efficiency and (| 6 l) 
assumes that energy harvesting from AWGN is negligible 
(very small compared to the desired signal ||9l)- In our analysis 
we assume C = 1 without loss of generality. In the second time 
slot, the received signal at the destination can be written as 


VD 




9r 


s -I- riD, 


(7) 


and the associated SNR is given by 


SNRd = 


Pr\9^? 

(1 -f c“)ct2 ■ 


( 8 ) 


In order to derive the outage probability of the RRS proto¬ 
col, firstly, we need to calculate the steady-state distribution 
of the battery. 

Proposition 1. The steady-state probability that a relay node 
is charged at the beginning of the broadcast phase, for the 
RRS scheme, is given by 




Jexp(-T') ^ 

5exp(-vl/)2i^ + ^- 


(9) 


Proof: See Appendix lAl ■ 

The expression in (|9]l shows that for highly dense networks 
A —oo, the steady-state probability that a relay is charged 
approaches to 1. This was expected, since as the density of the 
network increases, the relay selection probability approaches 
to zero ( 1 /A^ — 0 ) and therefore the relay nodes are most of 
the time in harvesting mode. On the other hand, the steady- 
state probability decreases as 'k increases is a decreasing 
function of T'); as the harvesting threshold increases (size of 
the battery), the probability to have an input power higher than 
the threshold decreases. The spectral efficiency does not affect 
the battery status of the relays, since uncharged relays observe 
only the energy content of the received signals, while charged 
relays are fully discharged in case of selection (independently 
of the decoding status (success or failure)). 

As for the outage probability performance of the RRS 
scheme, an outage event occurs when a) there are not any 
available relays in V i.e., N = 0, h) the disc contains at 
least one relay (N > 1) but the selected relay has an empty 
battery, c) TV > 1 , the selected relay is fully charged but cannot 
decode the source signal, and d) W > 1 , the selected relay is 
fully charged and decodes the source signal but the destination 
cannot support the targeted spectral efficiency. By analyzing 
the probability of these events, the outage probability for the 
RRS scheme is given by the following theorem. 

Theorem 1. The outage probability achieved by the RRS 
scheme is given by Coll. 

Proof: See Appendix IbI ■ 

Theorem [T] can be used to study the diversity gain of the 
RRS scheme, as shown in the following remark. 

Remark 1. For the special case with P —> oo, P^ —>■ oo, 
= Pr/P (constant ratio), p « do, and N > 1, the outage 
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IIrrs =exp(-A7rp^) + (l - exp(-A7rp^)) 


.-,r^l(^|£:)exp|-=|l + 


77 




exp ( — — + cIq — 2xdQ cos{9))i j xdxdd 




( 10 ) 


probability of the RRS scheme is given by the expression 


n 


oo 

RRS 


1 - s 



^ 1 - 


/ 1 + rfg 

I 



( ll a) 

( ll b) 


Proof: When N > 1, the first factor in (Uni) can be 
ignored, while Ci ~ do due to the assumption p « do. For 
P ^ oo (e.g., S —0), and by using the approximations 
1 — exp(—a;) r; x and 'y{a,x) ~ x°'/a for a; —>• 0 , the 
expression in Remark [1] can be obtained in a straightforward 
way. This remark assumes > 1 in order to highlight 
the impact of the battery on the achieved outage probability 
performance; it also holds for N >0 for moderate/high values 
of Xirp^. It is worth noting that the asymptotic expression in 
(fTTbl) is general and holds for any do. ■ 

The above expression shows that the outage probability for 
the RRS scheme suffers from an outage floor at high SNRs, 
which depends on the steady-state distribution of the battery. 


B. Relay selection based on the closest distance 

The relay selection based on the closest distance (RCS) 
requires an a priori knowledge of the location of the relay 
nodes. We assume that the AP monitors the location of 
the relays through a low-rate feedback channel or a Global 
Positioning System mechanism, and selects the relay node 
which is closest to the AP. The RCS scheme does not take 
into account battery status and/or instantaneous fading and also 
corresponds to a low implementation complexity, specifically 
for scenarios with low mobility. The relay node selected by 
the RCS policy can be expressed as 

R* = arg^.g,j, jnin^di. (12) 

The mathematical description of the RCS scheme follows the 
equations in Section IIII-AI by simply replacing the random 
relay node Ri with R*. The steady-state distribution of the 
battery for the RCS scheme is given by the following propo¬ 
sition. 


Proposition 2. The steady-state probability that a relay node 
is charged at the beginning of the broadcast phase, for the 
RCS protocol, is given by 


= vr- 


(13) 


Proof: Both RRS and RCS schemes select a single relay 
for transmission without taking into account the battery status. 
In addition, the selected relay is discharged at the end of the 
relaying slot independently of its decoding efficiency. Since 
both schemes handle the selected relay in the same way and 
the probability to select a relay in the RCS scheme is also 
E[l/A^] Ri I/Attp^ (i.e., a relay can be the closest with the 


same probability), the steady-state distribution of the RCS 
scheme is equivalent to the RRS scheme. ■ 

For the outage probability performance, the scenarios for 
an outage event follow the discussion in Section IIII-AI The 
following theorem provides an exact expression for the outage 
probability achieved by the RCS scheme. 

Theorem 2. The outage probability achieved by the RCS 
scheme is given by (dll. 

Proof: See Appendix ICl ■ 

In order to simplify the outage expression of the RCS 
scheme and study the diversity gain of the system, we provide 
the following remark. 

Remark 2. For the special case with P —^ oo, Pr —>■ oo, 
'^ = Pj-jP (constant ratio), p << do, N > I and a = 2, the 
outage probability of the RCS scheme is given by 


(15a) 

^ 1 - (15b) 

Proof: See Appendix iDl ■ 

The above remark shows that the outage probability of 
the RCS scheme converges to an outage floor at high SNRs, 
which depends on the charging behavior of the battery. By 
comparing the expressions in Remarks [T] and |2] we can see that 
both schemes converge to the same outage floor and therefore 
become asymptotically equivalent (the convergence floor is 
independent of a and do). However, by carefully comparing 
the expressions in (II lal) with a = 2 and (I15al) . it can be seen 
that the RCS scheme converges to the outage floor faster than 
the RRS scheme. 


U 


RCS 


1-pf 


RCS 


Att 


1 


Att 

„ f 1 -b 


- ( 1 + 


exp(—A ttp^) 

1 — exp(—A ttp^) 


-bl 


C. Random relay selection with battery information 

The random relay selection with battery information 
(RRSB) scheme randomly selects a relay node among the 
charged relays (if any). The RRSB scheme is based on a priori 
knowledge of the battery status and requires relays to feed 
their battery status ( 1 -bit feedback) at the beginning of each 
broadcast phas^ The steady-state distribution of the battery 
is given as follows 

Proposition 3. The steady-state probability that a relay node 
is charged at the beginning of the broadcast phase, for the 

^ Distributed implementation can be also considered where synchronized 
local timers at the relays (which could depend on the geographical location) 
allow the charged relays to access the channel according to the considered 
relay selection policy (28). 
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Hrcs = exp(-A7rp^) + “ exp(-A7rp^)^ 


Rcs 2 ’rA 2 exp(-H(l + i)) 


[1 — exp ( — 7rAp^)]^ 


r27v rp rf 

Jo Jo Jo 


exp f + ^0 ~ 2rdQ cos{9)) s — X-kt^ — — Xnx^j rxdrdxdO 


(14) 


=Q'(A) 


RRSB scheme, is given by 


= 1 - 


Vi/'S 

5\-k exp(—'I') 7 (^, '^p°‘) ’ 


(16) 


Proof: The charged relays form a PPP Q. which yields 
from the original PPP T* by applying a thinning operation; its 
density is equal to An = let N' denote the number 

of the charged relays in the disc, which follows a Poisson 
distribution with parameter An- The RRSB scheme randomly 
selects a single relay from the PPP fl and thus the probability 
that a fully charged battery becomes empty is equal to the 
selection probability i.e., 


^RRSB ^ ]g 


1 

W' 


1 

- E[N'] 


1 

An7rp2 


1 

A77®R®B7rp2 


(17) 


where the expression in ([TtT i applies Jensen’s inequality. On 
the other hand, the probability that an empty battery becomes 
fully charged (tto) is similar to the RRS scheme and follows 
the analysis in Appendix By substituting ttq and tt®®®® 
into (O and solving the linear equation for 77 ®®®®, we have 

77 ®R®b = with the required condition ttq > 1 /A 7 rp^. 

It is worth noting that the proposed approximation al¬ 
lows the derivation of a simple closed-form expression for 
the steady-state distribution and becomes more efficient for 
Xirp^ >> 0. The efficiency of the proposed approximation is 
discussed in Section IV] ■ 

By combining the expressions in (|9]l and (fTTl i. we can show 
that p®®® = 1/(2 — p®®®®) > ? 7 ®®®®. This observation is 
justified by the fact that the RRSB scheme limits the selection 
only between the charged relays and thus a charged relay can 
be discharged with a higher probability tti- 

As for the outage probability performance, an outage event 
occurs when a) there is not any charged relay in the disc, or 
b) the first or the second hop of the relay transmission is in 
outage. By analyzing these events, the following theorem is 
given. 


Theorem 3. The outage probability achieved by the RRSB 
scheme is given by the expression in 

nRRSB=exp(-An7rp^) + (1 - exp(-Ao 7 rp^))(l - Q), 

(18) 

where Q is the success probability for the relaying link and 
is defined in (doll. 


Proof: Based on the above outage events, the outage 
probability for the RRSB scheme is written as 

Hrrsb = P{iV' = 0} + P{Ar' > 1} 

X [1 - P{SNR, > €\N’ > l}P{SNRr, > e\N' > 1}] 
= exp(-An7rp^) + (1 - exp(-An7rp^))(l - Q)- 


When P,Pr 00 with = Pr/P (constant ratio), 
we can straightforwardly show that Q ^ 1 and thus the 
outage probability is dominated by the event where no relay 
is fully charged in the disc; the outage probability in (fTSl l 
asymptotically converges to the following outage floor 

n^SB exp(-An 7 rp 2 ) = exp(-Apf®®® 7 rp^). (19) 


D. Relay selection based on the closest distance with battery 
information 

The relay selection based on the closest distance with 
battery information (RCSB) scheme follows the principles of 
the RCS scheme, but takes into account the battery status of 
the relays nodes. Specifically, the RCSB scheme selects the 
closest charged relay according to 

min di. (20) 

As far as the steady-state probability of the battery is con¬ 
cerned, the charging/discharging behavior of the batteries 
follows the discussion of the RRSB scheme. 

Proposition 4. The steady-state probability that a relay node 
is charged at the beginning of the broadcast phase, for the 
RCSB scheme, is given by 

^RCSB^^RRSB ( 21 ) 

By straightforwardly applying the analysis of the RCS 
scheme for the outage events discussed in the RRSB scheme, 
the outage probability for the RCSB scheme is given by the 
following theorem. 

Theorem 4. The outage probability achieved by the RCSB 
scheme is given by the expression 

IIrcsb = exp(-An7rp^) + (1 - exp(-An7rp^))[l - Q'(An)], 

( 22 ) 

where Q'{-) is the success probability for the relaying link 
defined in dll. 

When P,Pr 00 with vk = P/Pr (constant ratio), by 
using the same arguments with dUl, we can show that the 
outage probability in (| 2 ^ asymptotically converges to the 
following outage floor 

n^SB exp(-An7rp^) = exp(-A77f'®®®7rp^). (23) 
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E. Distributed beamforming 

The distributed beamforming (DB) scheme selects all the 
relay nodes that are charged at the beginning of the broadcast 
phase; it is mainly used as a useful performance benchmark 
for the single-relay selection schemes. More specifically, all 
the relay nodes with fully charged batteries become active 
and attempt to decode the source signal. The relays that are 
able to successfully decode the source signal, create a virtual 
multiple antenna array and coherently transmit the signal to 
the destination (virtual multiple-input single-output channel). 
The practical implementation of the distributed beamforming 
requires the encoding of the source message with a cyclic 
redundancy check (CRC) code for error detection. In this way, 
only the relays whose CRCs check transmit in the second 
phase of the protocol. The DB scheme requires a perfect 
time synchronization and signaling between the relay nodes 
as well as a CSI at the relays; a feedback channel ensures the 
knowledge of the gi channel coefficient at the i-th relay. The 
broadcast phase of the DB protocol follows the description of 
the RRS and RCS schemes and is given by the expressions 
in (l4]i-(|6]l. The received signal at the destination, during the 
second phase of the protocol, can be written as 


VD 



WjQi 


s + no, 


(24) 


where Wi = g */precoding coefficient at 
the 7-th relay that ensures coherent combination of the relaying 
signals at the destination, and C denotes the set of the charged 
relay nodes (index), which successfully decoded the source 
message and participate in the relaying transmission (ED. The 
associated SNR at the destination is given by 


smD = PrY^ 

i&C 


(l-bc“)cr2- 


(25) 


Equivalently to the previous protocols, the outage probabil¬ 
ity performance of the DB scheme depends on the probability 
of a battery to be charged at the beginning of the broadcast 
phase. Towards this direction, we state the following proposi¬ 
tion. 


Proposition 5. The steady-state probability that a relay node 
is charged at the beginning of the broadcast phase, for the DB 
protocol, is given by 


S exp {-'if) ^ 
S exp (-T') 


(26) 


By plugging (l40l i and tti = 1 into rji of (O, Proposition |5] is 
proven. ■ 

The expression in (l26l l shows that the steady-state probabil¬ 
ity is independent of the density A. The DB scheme enforces 
all the charged relays to be fully discharged at the end of the 
relaying time slot and thus the size of the network does not 
affect the battery status distribution. 

Based on the description of the DB scheme, an outage 
event occurs when a) the relay set C is empty or b) when the 
coherent relaying transmission is in outage i.e., the destination 
is not able to decode the relaying signal. For the outage 
probability performance of the DB scheme, we state the 
following theorem. 

Theorem 5. The outage probability achieved by the DB 
scheme when p « do, is given by 

OO 

Hob = XI “ 

where A' = Ap™^ exp (-S) 

Proof: See Appendix [E] ■ 

It is worth noting that the sum in (l27l i quickly converges to 
the outage probability of the system and only a small number 
of terms is required (i.e., less than 10 terms). For the high 
SNR regime, we provide the following simplified expression, 
given as a remark. 

Remark 3. For the special case with P —^ oo, Pr —>■ oo, 
41 = Pr/P (constant ratio), the outage probability of the DB 
scheme is given by 

n“ « exp(-A'V)Jo 

—>■ exp(—Ap^Trp^). (28) 

Proof: See Appendix 10 ■ 

The expression in Remark |3] shows that for high SNRs, the 
outage probability of the DB scheme is equal to the probability 
that the set C is empty. This probability is an exponential 
function of p™ and approaches zero as Ap™p^ increases 
However, the DB scheme corresponds to a higher complexity, 
since it requires a continuous feedback to enable coherent 
combination of the relaying signals at the destination. 

Remark 4. For the special case with P —i oo, Pr —>■ oo, 
47 = Pr/P (constant ratio) and Xitp^ » 0, the outage 
probabilities of the proposed relay selection schemes are 
ordered as 




exp(—A' ttp^) 


2\k 


(A'V) 

k\ 


(27) 


Proof: The behavior of an empty battery follows the dis¬ 
cussion of the RRS and RCS schemes and therefore ttq is given 
by the expression in (l40l) . On the other hand, the proposed DB 
scheme enforces a fully charged battery to be discharged at 
the end of the relaying slot, since any charged relay is selected 
for potential transmission. According to the considered two- 
state battery model, any charged selected relay is discharged 
at the end of the relaying slot whether participating in the 
relaying transmission (successful decoding) or not. This means 
that the transition probability tti in this case is equal to one. 


nRcsB — n^RSB < — h^rs ■ (29) 

Proof: The proof follows straightforwardly from the 
expressions derived in (II Ibb . (I15bl) . ( fT9] l. (|2^ and (|28]) . 
Regarding the battery-based selection schemes, which result 
in an exponential outage probability, we have g™ = < 

^RRSB ^ I _ ^ ^ gjj.pj'_^j^RRSB^p2^ ^ exp(—Ary^TTp^) 

with 1 / = sexp{-^)-y{S ^ Att; due to path-loss attenua¬ 
tion, practical implementations consider small 4> and thus this 
condition is satisfied with strict inequality. The intuition behind 




















Fig. 3. Network topology for the multi-cell scenario; Voronoi cells are 
approximated by discs of radius l/4y^ and each relay transmits in a 
dedicated channel. 


this result is that in the asymptotic regime, where the achieved 
outage probability floor is a function of rji and independent of 
the channel, the selection of a single relay ensures connectivity 
while keeps more charged relays (and thus a higher rji) than 
the DB scheme. ■ 

IV. Generalization to multi-cell scenarios 

In this section, we discuss the generalization of the proposed 
relay selection schemes for large-scale networks with multiple 
APs. We focus on the single-relay selection schemes (RRS, 
RCS, RRSB, RCSB) and we study their behavior when multi¬ 
user interference affects the decoding operation at the relay 
nodes (well-known model in the literature e.g., ll^ , ISOl). 
More specifically, we assume a multi-cell network, where APs 
form a PPP T with density p on the plane K^. The relay nodes 
are also spatially distributed on the plane according to a 
PPP with density A. Each relay is connected to the closest AP 
and therefore all the relay nodes, which are located inside the 
same Voronoi cell, are dedicated to assist the corresponding 
AP. We study a typical cell, V, where the AP is located at the 
origin (Slyvnyak’s Theorem ll34ll ). For tractability, we assume 
that each Voronoi cell can be approximated by a disc of radius 
p = on and each AP has a single destination at 

Euclidean distance p in some random direction (worst case 
scenario). A direct link between AP and destination is not 
available and communication can be performed only through 
the DF relay nodes. In addition, we assume that each relay 
node transmits in an orthogonal channel and thus destinations 
are free from multi-user interference. This setup is inline with 
modern network architectures, where the relay nodes (e.g., 
femtocells, distributed antenna systems, etc.) have cognitive 
radio capabilities and thus can opportunistically access the 
channel. By using appropriate sensing radio mechanisms, the 
relay nodes exploit unoccupied spectrum holes or white spaces 
in order to minimize/mitigate interference ll^ . ll^ . The 
topology considered is an extension of the single-cell case and 
therefore the system model follows the discussion in Section 
Hn (except if some parameters are defined otherwise); Fig. [3 
depicts the network topology for the multi-cell case. 


A. RRS scheme 


According to the principles of the RRS scheme, each AP 
selects a relay node that is inside its Voronoi cell, in a random 
way. Since WPT highly depends on the Euclidean distance, 
we assume that WPT from external APs is negligible i.e., 
each relay node harvests energy only from the associated 
AP. Therefore, the multi-cell scenario does not change the 
analysis for the steady-state distribution of the battery and the 
probability that a relay node is charged, is given by Proposition 
[T] with p = 1 Let Ri denote the selected relay for the 

typical cell; the signal to interference and noise ratio (SINR) 
at the i-th relay is given by 


SINR, = 


P\hi\^ 

df 

pTT^’ 


(30) 


where / = J2jer/{v} 7^ denotes the aggregate (normalized) 
multi-user interference at the typical relay node, Hj denotes 
the channel power for the link between the j -th interfering AP 
and the selected relay, and rj denotes the associated Euclidean 
distance. It is worth noting, that for the broadcast phase of 
the cooperative protocol, the path-loss function 1 -f (i“ is 
replaced by the conventional unbounded model ll34ll : this 
assumption significantly simplifies the analysis and allows the 
derivation of closed form expressions. For the second phase 
of the cooperative protocol, since relay nodes transmit in 
orthogonal channels, the SNR at the destination is given by 
ill. For the outage probability of the RRS scheme, we state 
the following theorem. 


Theorem 6. The outage probability of the RRS scheme for the 
multi-cell case is given by 


IIrrso = A(A, 

(31) 

where A(-) is defined in (132b. 


Proof: See Appendix iGl 

■ 


For the high SNR regime, the outage probability of the RRS 
scheme can be simplified as follows 


Remark 5. For the special case with P —^ oo, Pr —>■ oo, 
Af = PrjP (constant ratio) and N >1, the outage probability 
of the RRS scheme is given by 

n^RSo^i-^rQi (0,1/4 Vm) (33) 

where Qi{-) is defined in (132b . 

In comparison to the single-cell case, we can see that the 
multi-user interference affects the first hop of the relaying 
protocol and therefore the achieved outage probability con¬ 
verges to a higher outage floor i.e., for the same radius p, 
n^RS = 1 - = 1 - ^rQi(0,p), where 

Qi(0,p) < 1 (see (l32]i) is the probability to successfully 
decode the source message at the relay node for high SNRs. 


B. RCS scheme 

The single-cell RCS scheme can be straightforwardly ex¬ 
tended to the multi-cell scenario by applying the RCS policy 
at each Voronoi cell. More specifically, each AP selects the 
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A(A, p) = exp(—A ttp^) + (l — exp(—A ttp^)) 



( 

119-0 1 l 


\ 

rexp(-Sx“)exp 

Jo 

TT 


ex^ 


16 


ex^ + p“ 

y 


xdx 


=Ql(H,p) 


TTp^ 


• “TiT 


r‘2'K rp 

J J exp -^(r^ + — 2rpcos(0))7 ] rdrd9 


(32) 


closest relay node among the relays which are located inside 
its circular coverage area. Equivalently to the RRS scheme, the 
multi-user interference affects only the broadcast phase of the 
cooperative protocol; therefore the steady state probability of 
the battery is equivalent to the single-cell case i.e., with 
p = 1/ 4^/p. For the outage probability, we state the following 
theorem. 

Theorem 7. The outage probability of the RCS scheme for 
the multi-cell case is given by 

HrcSo = 0(^1 1/4v/m), (34) 

where 0(-) is defined in (I351 l. 

Proof: The proof follows the analysis of the RCS scheme 
for the multi-cell case. It can be obtained in a straightforward 
way by using the PDF of the closest distance given in ( |46] | in 
order to calculate the probability to successfully decode the 
source message at the relay. ■ 

For the high SNR regime, the achieved outage probability 
can be simplihed as follows 

Remark 6. For the special case with P —> oo, Pr —> oo, 
'h = Pr/P (constant ratio) and N >1, the outage probability 
of the RCS scheme is given by 

nRcs„» = 1 - vT^Q'i (A, 0, l/4v/^) (36) 

where is defined in 

An interesting observation is that in contrast to the single¬ 
cell case, the RRS and RCS schemes do not converge to same 
outage floor at high SNRs. The RCS scheme converges to a 
lower outage floor, since the selection of the closest relay is 
a mechanism to protect the source signal against multi-user 
interference. 

C. RRSB/RCSB schemes 

Equivalently to the above discussion, the RRSB and RCSB 
schemes can be straightforwardly extended to the multi-cell 
scenario, by applying the RRSB and RCSB policies at each 
Voronoi cell. Since the generalization of the protocols does 
not modify the steady-state distribution of the battery and only 
affects the decoding probability at the relay nodes, the outage 
probability of the RRSB/RCSB is given by 

nRRSBo = A(Ap^®^l/4v/M), (37) 

IIrcsbo = 1/4^^). (38) 

By using similar arguments with the previous cases, 
the RRSB/RCSB schemes asymptotically converge to 


ng^RSBo = exp(-AprB) + (1 _ exp(-Apr®))(l - 
77r'^Qi(0,l/4VM)) and H^csbo = exp(-Apr«) + (1- 

exp(-A77fCSB))(l _ ,yRCSBQ/(^^RCSB^o^;^/4^))^ 

tively. 

V. Numerical results 

Fig. m plots the outage probability performance of the 
proposed relay selection schemes versus the transmitted power 
P. The first main observation is that the RRS and the RCS 
schemes converge to the same outage floor at high SNRs, as it 
has been reported in Remarks [1] and |2] respectively. However, 
the RCS scheme slightly outperforms the RRS scheme at the 
moderate SNRs and thus converges to the outage floor much 
faster. The RRSB and RCSB scheme, which take into account 
the battery status and avoid selection of uncharged relays, 
significantly improve the achieved performance and converge 
to the lowest outage floor; both schemes converge to the same 
outage floor at high SNRs. On the other hand, the DB scheme 
outperforms RRSB/RCSB schemes at low and moderate P. 
For these values, the transmission of the source signal by 
multiple relays through beamforming, boosts the SNR at the 
destination and improves the outage probability performance. 
In addition, it can be seen that as the density A increases, more 
relays participate in the relaying operation and therefore the 
gap between DB and single-relay selection schemes increases. 
In the same figure, we plot the theoretical derivations given in 
([Tol l. (O, (fTsT l. (l22li . and (l27l i. The expressions in (fTOl i. (fT4l) . 
([H, (El) refer to the exact performance and perfectly match 
with the simulation results. On the other hand, the expression 
in (El efficiently approximates the performance of the DB 
scheme; although (El holds for p « do, we can see that it 
is a tight approximation for the considered setup with do = 2p. 

Fig. H shows the impact of the system parameters do and 
4/ on the outage performance of the proposed relay selection 
schemes. As it can be seen in Fig. Ha), as the distance between 
the AP and the destination increases, the outage probability 
increases; a larger distance corresponds to a higher path-loss 
degradation. The DB scheme is more robust to the distance 
increase, since multiple relays cooperate through beamforming 
to overcome the path-loss attenuation. On the other hand. Fig. 
Hb) shows that the parameter 4/ is critical for the achieved 
performance of the system; this parameters characterizes the 
harvesting capability at the relay nodes as well as the available 
power for relaying. As it can be seen, a small 4* facilitates the 
harvesting operation since the harvesting threshold is low, but 
the available power for relaying is not able to ensure successful 
decoding at the destination. If the threshold 4^ is too high, the 
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©(A, p) = exp(—A ttp^) + ^1 — exp(—A ttp^)^ 




2tt\ 


1 — exp ( — 7rAp2) 


f 


exp 


- Xttx^ - 

16 


„ 2i^l 1,2 ;2-5 ; 
ex°‘ \ 1 +^ 




ex°^ + p“ 


xdx 


1 / 


A exp -^ 


r2'K rp 

Jo Jo 


1 — exp ( — wXp^') 


^Q[(X,S,p) 


exp ( ~ 2rpcos(6)) « — X-wr^ j rdrdB 


(35) 



Fig. 4. Outage probability versus P\ p = 3m, do = 2p, 'I' = 0.1, a = 3, 
(T^ = 1, ro = 0.01 BPCU, and A = {1, 0.5}; the dashed lines represent the 
theoretical results. 


relay nodes transmit with a high power but the probability to 
be charged is decreased. 

In Fig. |6] we plot the outage probability performance 
for different spectral efficiencies rg- As it can be seen, the 
convergence outage floor of the relay selection schemes is 
independent of the spectral efficiency. Increasing the spectral 
efficiency affects only the convergence rate of the selection 
schemes (slower convergence). This observation has been 
expected, since according to our analysis, the convergence 
floor only depends on the steady-state distribution of the 
battery, which is not a function of the spectral efficiency. For 
the sake of presentation, we use small values of rg in our 
simulation results without loss of generality. 

Fig. |2] shows the impact of the radius p on the achieved 
outage probability performance. The first main observation is 
that as p increases, the performance of the proposed relay 
selection schemes is improved at high SNRs i.e., they converge 
to lower outage floors. This observation was expected, since 
as the area of the disc increases, the number of the relay nodes 
increases (with E[A^] = Xirp^) and thus a) more relays partici¬ 
pate in the relaying operation (DB scheme), b) the probability 
to select an uncharged relay node is decreased (RRS, RCS 
schemes), c) the probability to have no charged relay in the 
disc is decreased (RRSB, RCSB schemes). However, for the 



Fig. 5. Outage probability versus a) do and b) Simulation parameters: 
P = 30 dB, p = 3 m, a = 3, 0-2 = 1, ro = 0.01 BPCU, A = 0.5, a) 
'P = 0.1, b) do = 2p. 



P[dB] 

Fig. 6. Outage probability versus P for different values of ro- Simulation 
parameters do = p, '^ = 0.1, a = 3, = 1, ro = {0.1,1.5} BPCU, 

A = 1; the dashed lines represent the theoretical results. 
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Fig. 7. Outage probability versus P for different values of p. Simulation 
parameters: do = p, ^ = 0.1, a = 3, <j^ = 1, vq = 0.1 BPCU, A = 1; the 
dashed lines represent the theoretical results. 



Fig. 8. Steady state-distribution of the battery for the RRSB/RCSB schemes; 
efficiency of the proposed approximation in GD for different settings (A, p). 


low and moderate SNRs, it can be seen that as p increases, the 
outage probability performance decreases. In this SNR regime, 
the outage probability is dominated by the channel path-loss 
attenuation and becomes more severe as p increases; a larger 
p corresponds to a higher distance (do = p) between source- 
destination and therefore increases the longest hop link. 

Fig. [8] deals with the steady-state distribution of the battery 
in the RRSB and RCSB schemes. The closed form expression 
given in Proposition |3] is based on the transition probability 
.^RRSB, yvjjjQh ^sgs fjjg Jensen’s approximation (inequality); see 
the proof of Proposition [3 In order to show the efficiency of 
this approximation. Fig. [8] compares the exact (given 

by simulation results) against the proposed approximation for 
different system parameters. As it can be seen, the proposed 
expression efficiently approximates the steady-state distribu¬ 



Fig. 9. Outage probability for the RRSB and the batteryless RRS 1181 
schemes versus P. Simulation parameters: do = p, = 0.02, o: = 3, 
0-2 = 1 = 0.01 BPCU, A = 1. 



Fig. 10. Outage probability versus P. Simulation parameters: do = p, ^ = 
0.1, Q = 3, 0-2 = 1, ro = 0.001 BPCU, A = 0.5, and /j. = {0.001, 0.005}; 
the dashed lines represent the theoretical results. 

tion and its accuracy is improved as the term ttq — 1/A7rp^ 
increases. Therefore, this approximation allows closed-form 
expression for the steady-state distribution and provides accu¬ 
rate results for the scenarios of interest i.e., Xnp^ >> 0. 

In Fig.|9l we compare the proposed RRSB scheme with the 
battery less RRS scheme proposed in ifTSll . The battery less RRS 
scheme uses the PS-SWIPT technique at the relay node and the 
harvested energy is directly used to power the relaying link. 
It can be seen, that the RRSB scheme outperforms batteryless 
RRS at low and moderate SNRs and the gain increases as 
the conversion efficiency ( decreases and/or the disc radius 
increases. For this SNR regime (which is expanded as de¬ 
creases), the PS-SWIPT technique becomes inefficient and the 
integration of the battery, which decouples the information and 
power transfer in time, significantly improves the performance. 
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For high SNRs, the RRSB scheme suffers from an outage 
floor, while batteryless RRS provides a diversity gain equal to 
one ifTSll . It is worth noting that the batteryless RRS scheme 
assumes perfect channel knowledge at the relay node, and 
requires appropriate electronic circuits in order to perform 
dynamic PS; therefore it corresponds to a higher complexity 
than the proposed RRSB scheme. 

Fig. [TOl deals with the application of the proposed single¬ 
relay selection schemes to a multi-cell scenario with multi-user 
interference at the relay nodes (Section HVl i. The RCS scheme 
outperforms the RRS scheme for the low/intermediate SNRs, 
while both schemes converge to similar outage probability 
floors; for this simulation setup, the success probabilities at 
the relay nodes become almost equal at high SNRs for both 
p, values. The consideration of the battery status into relay 
selection significantly improves the outage performance of the 
system and thus RRSB and RCSB schemes achieve lower 
outage probability floors. It is worth noting that RRSB and 
RCSB schemes do not converge to the same outage floor at 
high SNRs, since the decoding probability at the relays is 
different for the two schemes and affects their convergence 
(see Section ITV-Cl l. On the other hand, as the AP’s density p 
increases, the radius of the Voronoi cells decrease and the 
curves follow the discussion in Fig. [T] Theoretical results 
perfectly match with the simulation curves and validate our 
analysis. 


VI. Conclusion 

In this paper, we have studied the problem of relay selection 
in WPC cooperative networks with spatially random relays. 
We assume that the relay nodes are equipped with batteries 
and use the received signal either for conventional decoding 
or battery charging. Based on a single-cell network topology, 
we investigate several relay selection schemes with different 
complexities. Their outage probability performance is derived 
in closed form by modeling the behavior of the battery as 
a two-state MC. We prove that the relay selection schemes 
suffer from an outage floor at high SNRs, which highly 
depends on the steady-state distribution of the battery. The 
RRS and RCS schemes achieve the worst outage probability 
performance and converge to the same outage floor, while 
the consideration of the battery status significantly improves 
their achieved performance. The DB scheme is a promising 
solution for low/moderate SNRs and outperforms single-relay 
selection at the cost of a CSI. The proposed selection schemes 
are generalized to multi-cell network topologies, where multi¬ 
user interference affects relay decoding. 


Ffi{x) = 1 — exp(—x). The CDF of the random variable Ui is 
given as follows 


Fu{x) = P{uj < a;} 

= l-P{|/i,P >x(l + dn} 

= 1 — Eexp {—a:(l -f df)} 

= 1-/" exp{-x(l-f y“)}/d(y)dy (39a) 

J-D 


= 1 - 


IV 

1 


f'2‘K np 


Jo Jo 


exp {—a;(l -I- y°‘)} ydydO 


2 

= 1-j exp {-x) / y exp {-xy°‘) dy 

P Jo 

5 -f{S,xp°‘) 

= 1-2 (“2;)-j-, 

x° 


(39b) 


where fd{x) = in ( I39ab denotes the probability density 

function (PDF) of each point in the disk T) and the result in 
(I39bl i is based on 1261 Eq. 3.381.8]. 

By using the CDF of the random variable Ui in ( l39l l. the 
probability that an empty battery is fully charged during the 
broadcast phase can be expressed as 


TTO 



— exp(-T') 
P 




(40) 


On the other hand, if N is the number of relays in D, the 
probability to select a relay according to the RRS policy is 
equal to 1/W; the probability that a charged relay becomes 
uncharged is equal to this selection probability and can be 
expressed as 


RRS 


= E 




= exp(-A7rp^)^ 


2\k 




(A^) 

fc! 


> 


E[N] A7rp2 


(41) 


where (HTt is based on Jensen’s inequality and E[A^] = Avrp^ 
is the average number of relays in V (from the definition of a 
PPP). It is worth noting that the probability is simplified 
by using Jensen’s inequality; this approximation significantly 
simplifies our derivations and is tight as Avrp^ increases. By 
combining (l40t and dril l and substituting back into dJll, the 
expression in Proposition [1] is proven. 


B. Proof of Theorem |7] 

The outage probability for the RRS scheme can be ex¬ 
pressed as 


Appendix 

A. Proof of Proposition [7] 

Firstly, we calculate the cumulative distribution function 
(CDF) of the random variable Ui = this result is 

essential for the derivation of the steady-state distribution. 
From the system model, \hi\'^ is an exponential random 
variable with unit variance and thus its CDF is equal to 


Drrs = P{iV = 0} + P{iV > 1, S{R,) = so} 

-b ¥{N > 1, S{Ri) = si, SNR, < ej 

-b F{N > 1, S{Ri) = si, SNR, > e, SNR^, < ej 

= F{N = 0} -b F{N > 1}(1 - + P{W > 


X 1 - P{SNR, > e|At > l}P{SNRu > e|At > 1} 

Q 

= P{W = 0} -b F{N > 1}(1 - 


(42) 













13 


From the Poisson distribution, we have 


at the destination, conditioned on > 1, is expressed as 


F{N = 0} = exp(-A7rp2). (43) 

In addition, 

P{SNR, > e\N > 1} = F{u, > S} = 1 - F,(S), (44) 

where Fu{-) is given in (|39] |. Finally, the success probability 
for the relaying link is written as 


{SNRij > e|7V> 1} =P<{ > 

(1 + cf.)' 


4- 


= E 


jexp 


(48a) 


41 

A exp (-S/^i) 

1 — exp ( — Attp^) 

2tt i-p / '^ 

expl —+ dp — 2rdo cos{0))^ — Xirr'^ 1 rdrdO. 


0 JO 


4- 


(48b) 


P{SNRi 5 > e|A^ > 1} = P \gif > 


-(1 + c “) 


= E exp ( — 


^(1 + cf) 


= E exp ( ~ ^ + (di + dQ — 2dido cos{6 ))« 


(45a) 


= / exp - — 


Jv 

1 

TTp^ 

X / 


41 L 


■exp - — 


1 + + dg — 2xdo cos{0)) j fdi{x)dx 


r^TT j-p 

0 Jo 


4- 


exp ( —+do — 2xdocos(0))^ ) xdxdO, 


(45b) 


where cf = df + dg — 2dido cos(d) in (145al) holds by using 
the cosine law. By combining (021, (04l) . (02 and substituting 
back into (02- we prove Theorem [1] 


C. Proof of Theorem |2 

The PDF of the nearest distance di- for the homogeneous 
PPP $ with intensity A, conditioned on TV > 1, is given by 
lE] Eq. (33)] 

2A7r 

fr{r) = - - --—— r exp ( - A7rr^). (46) 

1 — exp ( — Airp^j 


The outage probability for the RCS scheme can be ex¬ 
pressed by the general expression in (02l l. By combining (|2, 
(03l l. (07l i. (I48bb and substituting into (02l l. we prove the 
statement in Theorem [2] 


D. Proof of Remark^ 

For the special case that a = 2, ( 071 ) is simplified to 


‘{SNR, > e|iV > 1} 
27rAexp(—5) f 
1 - exp ( - 7rAp2{ Jg 
Att exp (—5) 

1 — exp ( — 7rAp^) 
A7r(l -S) 


Att + [ 




^ 1 + Sp 


rexp (—(■=. + A7r)r^) dr 

1 - exp ( - (S + \tt)p^) 
(5 + Att) 
exp(—A ttp^) 


1 — exp(—A7rp2) 


(for : 


(49a) 

^0), 

(49b) 


where ( 149 ab holds from Eq. 3.381.8], the asymptotic 
expression in ( I49bb uses the approximation 1 — exp(—x) ~ x 
for X —>■ 0. 

As for the relaying link, for the special case that p << do, 
the Euclidean distance of the link relay-destination becomes 
a* Ri dg and the expression in (148 ab is simplified to 


'{SNRd > e|A^> 1} = exp | — 
r; 1 - 


:(1 -f dg) 


(1 -f dg) 




(for 


0 ). 

(50) 


The probability to successfully decode the source message By combining (I49bb . (l50b and substituting into (02 with 
at the selected relay, conditioned on > 1, can be expressed P{A^ = 0} = 0, we prove Remark |2 
as follows 


P (SNR, > e|TV > 1} = P > e| 

= P{|/r,.|^>S(l + d“.)} 


= E{exp(-^(l-f d“.))} 

= [ exp(-5(l-f r“))/,.(r)dr 
Jo 

27rAexp (—S) 

1 — exp ( — TrXp'^) 


exp ( — Sr“ — Xirr^ I rdr. (47) 


The probability to successfully decode the relaying signal 


E. Proof of Theorem 0] 

In the DB scheme, the relay nodes which are active (charged 
batteries) and are able to decode the source message, par¬ 
ticipate in the relaying transmission. The participating relays 
{i € C) form a PPP <!>', which yields from the original 
PPP $ by applying an independent thinning operation. More 
specifically, the PPP has a density 

A' = Ap“(l-F„(5)). (51) 

In order to simplify the analysis, we focus on a disk V with 
p « dg; in this case, the distance between each relay and 
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the destination becomes equal to do i.e., Ci « dp. The SNR 
expression in ( |25] | is simplified to 

The random variable Y = X^ieC of 1^1 = 

K independent and identically distributed (i.i.d.) exponential 
random variables, where the cardinality K follows a Poisson 
distribution with density A'. Therefore, the outage probability 
of the system is given by: 


IIdb = PjSNR^) < e} 

■(1+d^) 


= p<^ r < 




= ^P{F < zo|iT = fe}P{iT = k} 


fc=0 


= E 


lik.zo) _^ 


A;=0 


nk) 


■ exp(—AVp^)- 


k\ 


(53) 


where FyiXyk) = denotes the CDF of the random 

variable Y i.e.. Gamma distribution with shape parameter 
k. By plugging (ISTl i into (l53T l. we prove the expression in 
Theorem |5] 


F. Proof of Remark]^ 

In the case that zq 0 (i.e., high SNRs with P oo), the 
expression in (l53] l is simplified to 


n 


OO 

DB 


exp(-A'7rp^)^ 

/c=0 


(zo)'=(A'7rp2)fe 

(fe !)2 


= exp(-A' 7 rp^)Jo ^ 2 p\/ zqA'tt^ 
exp(-Ap?®7rp2), 


(54a) 

(54b) 

(54c) 


where /o(-) denotes the zeroth order Modified Bessel function 
of the first kind, (I54ab uses the approximation 7 ( 0 , x) ~ x°‘/a 
for small x, (I54bb is based on the series representation in ll^ 
Eq. 8.447.1], and (I54cb uses the approximation Io{x) ~ 1 for 
small X. 


G. Proof of Theorem |6] 

The outage probability of the RRS scheme follows the 
analysis of the single-cell case, presented in Appendix iBl The 
main difference is the computation of the decoding probability 
at the relay node. This probability conditioned on di (distance 
between AP and relay), can be expressed as follows 

r p\hd^ 

P{SINR, > e|iV > 1, dj = P i „ > e 

PI + 

= P{|/i,|2 > 5d“ + e(i“/} 

= exp i-Ed“) exp {-ed’^I) 

= exp (-Sd^ L(ed“), (55) 


where L(-) is the Laplace transform of the random variable 
I = /{V} With expectation over di, we obtain 


n p 

exp{-Ex°‘)L{ex°‘)fd{x)xdxd9 

2 

= —^ / exp (—Sa;“) J 7 (ea;“)a;da;, 

P Jo 


(56) 


where fd{x) = denotes the PDF of each point in the 

(typical) disc. 

The Laplace transform of the interference I can be calcu¬ 
lated as follows 


L{s) = Eexp(—sJ) 


= E exp(-siJjr^. “) 

UGT/{D} 


= exp < — 27 r/j, / Ei/ 


1 — exp(—sJJr “) 



(57) 

where ( l57b uses the probability generating functional of a PPP 
im Sec. 4.6]. Conditioned on H, we have 


nOO 

(1 - exp(-siJr"“))rdr 
J O 


= 5 

= 5 


= exp - 


1 - exp 

1 — exp{—sHx)^x~^~^dx 
sH 


(58a) 

(58b) 


— 1 ) -I- sH J x ° exp{—sHx)dx 

(58c) 


= exp - 


sH 


-1] + {sHYj 1 - d, 


sH 


(58d) 


where ( I58ab follows from the substitution y r°‘, ( I58bb from 
the substitution x 4 — y~^, ( I58cb from integration by parts, and 
( I58db is based on ||26] Eq. 3.381.8]. With the expectation over 
H, we have 

exp(-JJ)dJJ 


P exp 


+ s" 


jS I - . S/i 


T ( 1 - exp(-JJ) ) dJJ 


1-5 2 F 1 


(1,2 ; 


2-d; 






s + p° 


(59) 

where ( l59b is based on the expression in ll2^ Eq. 6.455.2]. By 
substituting (l59b into (l57t . we have 

/ 


i(s)=exp 


-Tip 


V 


Fi(l,2;2-d;. 


r + 1 


r,a-2 




2 

p S 

s + p° 


(60) 
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By combining (|9]l, (l4?b . (l45T l, (l56T l and substituting back into 
(|42]i with do = p = 1/ 4:y/JI we prove Theorem |6] 
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